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ABSTRACT 

We present Atacama Large Millimeter/submillimeter Array observations of 99.02 GHz free- 
free and H40 g emission from the centre of the nearby starburst galaxy NGC 253. We calculate 
electron temperatures of 3700-4500 K for the photoionized gas, which agrees with previous 
measurements. We measure a photoionizing photon production rate of(3.2 ± 0.2) x 10^^ s“^ 
and a star formation rate of 1.73 ± 0.12 Mq yr“^ within the central 20x 10 arcsec, which 
fall within the broad range of measurements from previous millimetre and radio observations 
but which are better constrained. We also demonstrate that the dust opacities are ^3 dex 
higher than inferred from previous near-infrared data, which illustrates the benefits of using 
millimetre star formation tracers in very dusty sources. 
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1 INTRODUCTION 

The millimetre waveband contains two star formation tracers that 
directly trace photoionizing light from star forming regions while 
not suffering from dust obscuration effects, making it a supe- 
rior alternative to st ar fo rmation tracers in oth er wavebands (see 
ICalzetti et al.l (l2009h and iMurphv et al.l (l201lh for reviews). The 
continuum emission includes free-free emission, and while thermal 
dust emission is dominant at >100 GHz and synchrotron emission 
is more prominent at <30 GHz, free-free emission is the domi¬ 
nant emission source at 30-100 GHz (e.g jPeel et al.llimh . Mul- 
tiple hydrogen recombination lines are seen at millimetre wave¬ 
lengths as well. Optical and near-infrared lines are affected by dust 
extinction effects, and centimetre and metre recombination lines 
are affected by potential masing and collisional broadening effects, 
but millimetre lines ye not affected by any of these problems (e.g. 
iGordon & Walmsl^ll990h and are therefore more accurate tracers 
of star formation. 

The Atacama Large Millimeter/submillimeter Array (ALMA) 
is a groundbreaking millimetre and submillimetre telescope with 
the capability of reaching sensitivity levels in the 30-950 GHz 
rang e at least an ord er of magnitude better than other telescopes 
(see iLundgrenI l2013l for a technical overview). ALMA is capa¬ 
ble of detect free-free and recombi nation line emission from man y 
nearby infrared-luminous galaxies (IScoville & Murchikovall2013h . 
It is therefore of interest to use early ALMA observations to explore 
the telescope’s capabilities to detect these star formation tracers. 

In this paper, we present ALMA observations of free-free 
continuum emission and H40a emission at 99.02 GHz from the 
centre of NGC 253 that we use to derive star formation rates 
(SFRs) as well as the electron temperature (Te) of the gas. In 


additi on, we use the H 4 0a da ta together with near-infrared data 
from Engelbracht et al. (1 1998 ) to exam ine t he dust attenuation . 
iBolatto et" (120131) 7 Leroy et al.1 (l2015h . and iMeier et al.l (l2015h 
previously published analyses based on these data, but except for 
a brief mention of recombination lines in the latter paper, they fo¬ 
cused on the molecular gas, whereas we will concentrate on the 
photoionized gas. 

Radio and millimetre recombination lines have pre v iously 
been detected from NGC 253 (e.g. Seaquist & Belli \ 911 _\ 
Mebold etaP Il980l; lAnanthyamaiah & GossI 1 19961: IPuxley et alJ 


_ _ _ _ __:ley r,, , 

19971: iMohan et all l200l l2005l: iRodrieuez-Rico et all l2006l: 
Keole^et^ 201 ih. Th is includes the detection of H40a emission 


Puxlev et al. ( 1997h . However, our analysis will focus on mil- 


by 

limetre lines that are less susceptible to masing effects than most 
previous recombination line observations, which have been at lower 
frequencies. Additionally, most of the continuum at the location of 
the H40a line is free-free emission, allowing us to use the line 
and continuum data together to measure Te, whereas the lower 
frequency data are more affected by synchrotron emission. The 
ALMA data also have superior spatial resolutions, spectral reso¬ 
lutions, and sensitivities than most previously-published NGC 253 
data. This ultimately serves as an example of the future work on ex- 
tragalactic star formation in nearby galaxies that can be done with 
ALMA. 


2 DATA 

The data were originally acquire d as part of _ program 

2011.0.00172.S and published by (iBolatto et al.l l2013h . The 
observations covering the H40a line are from three execution 
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Figure 1. Images of the continuum surface brightness, the H40a intensity, and the H40 q; mean velocity in the central 32x32 arcsec of NGC 253. The velocity, 
which is relative to the Solar System Barycentre, is only shown for data where the H40a intensity is detected at the bcr level. The while oval at the bottom 
right of each panel shows the FWHM of the beam. The green regions in the continuum image show the total (T), east (E), central (C), and west (W) regions 
within which fluxes and spectra were measured. The spectra are shown in Figure[2| and measured quantities are listed in Table[T] 



blocks (EBs) performed on 07 May, 01 July, and 02 July 2012. 
Three locations along the major axis spaced by 25 arcsec were 
observed, and the total integration time per location was 29 min. 
Only 14-20 antennas were operational. The spectral window 
covering the H40a line was centered at a sky frequency of 
98.54 GHz, contained 3840 channels each with a width of 488 kHz 
(1.5 km s“^), and included both polarisations. Uranus was the flux 
calibrator, J2333-237 was the bandpass calibrator, and J0137-245 
was the phase calibrator. 

We reprocessed the data using the Common Astronomy Soft¬ 
ware Applications version 4.2.2, which includes a version of the 
flux calibration (Butler-JPL-Horizons 2012) that is more up-to-date 
than the one used to create the data in the ALMA archive. The 
visibility data were processed through steps to flag bad data and 
to calibrate the phase and amplitude as a function of frequency 
and time. Next, we rescaled the amplitudes of the visibility data 
so that the signal measured for JO137-245 in the individual execu¬ 
tion blocks matched the weighted average of these values, and then 
we concatenated the data. We produced two versions of the image 
cube (with and without continuum) using a clean algorithm with 
natural weightings in an interactive mode. Although the primary 
beam is 63 arcsec, we created image cubes covering a 32x32 arc¬ 
sec region (with 0.1 arcsec pixels) because no detectable emission 
is found outside this region. The cubes used in the analysis cover 
a sky frequency range between 98.80 and 99.09 GHz, and each 
channel in each cube covers 8 channels in the visibility data, which 
corresponds to 3.9 MHz (~11.8 km s“^). For display purposes, we 
also produced a cube that covered 97.64-99.44 GHz, which was the 
full usable range of the spectral window containing the H40a line, 
with the same spatial and spectral resolution. The full-width at half 
maximum (FWHM) of the reconstructed beam is 1.9 x 1.6 arc¬ 
sec; assuming a distance of 3. 44 zb 0.26 Mpc (based on the av¬ 
erage of distances measured bv IPalcanton et al.ll2009h . this corre¬ 
sponds to spatial scales of ~ 30 pc. The flux calibrat ion uncer¬ 
tainty reported in the ALMA Technical Handbook is 5% (iLundgrenI 
I 2 OI 3 }) . The 98.54 GHz flux densities for the bandpass calibrator 
(1.15zb0.04Jy) and the phase calibrator (1.09zb0.06Jy) have rela¬ 
tive uncertainties (based on the variations in measurements from in¬ 
dividual EBs) that are consistent with the Technical Handbook, al¬ 
though unidentifled systematic amplitude calibration effects could 
affect the flux densities. Cleaned image cubes of the bandpass and 
phase calibrators show <1% variations in the flux density from 
channel to channel, and no spectral features are seen at the position 


of the H40a line in the visibility data for the calibration sources, 
which indicates that any line emission > 1% of the continuum emis¬ 
sion is not caused by calibration issues. 

We had access to data from program 2011.0.00061.S (PI: 
Takano) that also covers the H40a line. Because the spectral set¬ 
tings were different from program 2011.0.00172.S, we did not use 
the data in our analysis, but we did examine the processed archival 
data for program 201L0.0006LS to confirm that the H40a emis¬ 
sion is detected in those data. 

H40a fluxes, mean velocities (relative to the Solar System 
Barycentre, and velocity FWHM were measured in the continuum- 
subtracted image cube, and the continuum at 99.02 GHz was mea¬ 
sured in the other image cube by fitting and removing the line 
emission. Images of the the continuum emission (based on all con¬ 
tinuum data within the spectral window), the H40a flux, and the 
H40a mean velocity are shown in Figure [T] Spectra integrated over 
the central 20 x 10 arcsec as well as integrated within regions cover¬ 
ing the three brightest regions are shown in Figure [2| Details of the 
measurements in these regions are listed in Table [T] The uncertain¬ 
ties in the measurements incorporate the noise per channel values 
listed in Table [TJ but except for the east (E) region, the accuracy 
in the continuum and H40a fluxes is primarily limited by the cali¬ 
bration uncertainties. Most of the signal outside the 20 x 10 arcsec 
ellipse is detected at <3cr, and we lack the uv coverage to recover 
>20 arcsec structures. 

The central star forming region can be divided into three knots 
that lie along the major axis of the galaxy. The central (C) and 
west (W) sources are detected at >10cr in both the continuum 
and H40a emission. The eastern source is detected at >10cr in 
continuum emission as well, but the peak is only detected at the 
~8cr level in H40a emission. Although all thre e sources had previ¬ 
ous been detected i n radio continuum emission (Mohan et al . 12005 I : 
iKenlev et alJl201 ih . only jAnantharamaiah & GossI 1996h had pre- 
viously shown radio recombination line emission from all three 
sources. The sam e struc ture is also seen in the Br 7 image from 
lEngelbracht et akl (Il998h . but the source is only marginally re¬ 
solved and does not appear to show the same peak in recombination 
line emission. 

The H40a emission from the central source is relatively 
broad (with a FWHM of 191 ± 4 km s“^) and appears to 
be rotating orthogonally to the plane of the galaxy, which 
is co nsistent with previous radio recombination line observa- 
tions (lAnantharamaiah & GossI 1 199^ : iRodrfguez-Rico et al.ll200d : 
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iKenlev et alJl201 ih . This may indicate that material has recently 
fallen into the central region, that the central source comprises two 
superimposed photoionized regions, or that the ionized gas is in 
a superwind that is flowing in a direction aligned with the minor 
axis. It is unlikely that an active galactic nucleus is present, as the 
H40a line would appear >1000 km s in width. We do not de¬ 
tect any asymmetry in the line emission from the central source as 
had been reported i n some prior recom bination line observations 
(iMohan et al.ll2002l : IKcdIcv et al.ll201lh but not others. The reason 
for this discrepancy is unclear, although one possibility is that the 
line shapes differ because the lower-frequency lines are probing 
less dense gas than the H40a line. 


3 ANALYSIS 


For determining the relative contribution of free-free emission to 
the 99.02 GHz continuum, the frequency coverage of the ALMA 
data by its elf is insufflcient for modelling the spectral energy 
distribution. IPeel et al] (l201lh indicate that most of the globally- 
integrated 99.02 GHz emis sion would be continuum. However, 
iRodrfguez-Rico et III (l2006h found a much smaller contribution of 
free-free emission to the central 30 arcsec, but they assumed that 
the slope of the synchrotron emission is flxed by the 5 and 15 GHz 
data points. If we instead repeat this fit but allow the slope of the 
synchrotron component to vary (as shown in Figure 0 we And that 
70 zb 10% of the 99.02 GHz emission is from free-free emission. 
We applied this correction to the measured continuum emission in 
the rest of our analysis. 

The ratio of the H40a flux density integrated over 

velocity v to the free-free flux density [cont^ can be written as 


f fj^(line)dv f Jy 
fu{cont) [jykms 


5.06 X lO*" 


€1/ 

\ 1 

-0.83 

'Te' 

erg s“^ cm~^ 

L GHzJ 


K 


( 1 ) 


based on equations from IScoville & Murchikoval (l2013h . The emis- 
sivity Cjy varies slowly as a function of electron density rie but 
varies strongly with Te. As shown in Table [T] we measure line-to- 
continuum ratios of 42.5 in the central region and ~50 in the east 
and west regions. We used these ratios and interpolated between 
the Ciy values from IStorev & Hummej (Il995h for case B recombi¬ 
nation and ne=10^cm“^to calculate Te for each of the regions in 
our analysis. The Te values of 3700-4500 K we measure are slightly 
low er than the previously-reporte d values from dPuxlev et al.l 19971) 
and iRodrfguez-Rico et al.l (l2006l) . but only by <lcr. The Te values 
are also comparable wit h values measured by IShaver et~aD (Il983h 
and IPaladini et al.l(l2004l) within the central 8 kpc of the Milky Way. 
If we did not apply the correction of 0.70 to the continuum emis¬ 
sion, the measured Te would increase to values of 5800-7100 K, 
which are still consistent with Milky Way values. If 50% of the 
continuum emission is free-free, then the resulting Te fall within 
the implausible range of 2500-3000 K. 

The SFR that we calculate depends upon both Te and the star 
formation hist ory. The presence of both photoionized gas and su¬ 
pernovae fe.g. lLenc & Tingavll2006l : IPumpadarath et al.ll20i3) im¬ 
ply that the star formation is continuous and has been ongoing 
for >5 Myr. We used v ersion 7.0.1 of the Starburst99 mod¬ 
els (iLeitherer et al.lll9^ to determine how to convert the number 
of ionizing photons produced as a function of time (Q) to SFR. 
This version of Starburst99 incorporates versions of the Geneva 
stellar evolution tracks that incorporate rotation, with the rotation 



98 


99 
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Figure 2. Plots of the spectra measured in the 20x 10 arcsec elliptical region 
(labelled as “Total”) and in the east, central, and west subregions identifled 
in Figure [T] The frequencies are rest frequencies based on the velocities in 
Table[T] The red lines show major spectral features identifled using Splat- 
ALOGUE^. The blue bar shows the noise per channel as listed in Table[T] 


velocities of zero rotation and 4 0% of the break-up v elocity rep¬ 
resenting two potential extremes (ILeitherer et al.lbofj) . Using the 
average of the results from these two versions of the Geneva tracks 
for Z=0.040 metal licity (the high est metallicity for which tracks 
are available) and a lKroupal (l2002h initial mass function for a mass 
range of 0.1-100 M©, we And that SFR = 1 M© yr“^ corresponds 
to Q = 1.85 X 10^^ s“^ after ~5 Myr of continuous star formatior0. 

Given this, SFR can be calculated from either the free-free or 


^ http ://www. cv. nrao. edu/php/splat/ 

^ The Q from the two rotation scenarios differ by a factor of ~ 1.8. Using 
Z=0.014 metallicity (the next lower metallicity for which Geneva tracks are 
available) had a <5% effect on the result s. Using o l der st ellar population 
models that did not include stellar rotation. [Kennicud Il998l) r eported a con¬ 
versio n factor that would increase SFR by a factor of ~2, and lMumhv et all 
(ioii gave a conversion factor that would increase SFR by 35%. 
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Table 1. Measurements of the continuum and H40a emission. 


Region 

Noise 

99.02 GHz 

H40o 

H40o 

H40o 

H40o 

Electron 

SFR 



per 

Free-Free Flux 

Flux 

Mean 

FWHM 

/Free-Free 

Temperature^ 

Free-Free^ 

H40a 


Channel® 

Density^ 


Velocity 


Ratio^ 





(mJy) 

(mJy) 

(Jykms-i) 

(km s“^) 

(km s“^) 

(km s“^) 

(K) 

(M© yr-l) 

(M© yr-l) 

Total 

5.2 

146.5 ± 17.2 

6.91 ±0.80 

244 ±2 

169 ± 10 

47.2 ±2.4 

3900 ± 300 

1.59 ±0.16 

1.87 ±0.18 

East 

1.3 

13.4 ± 1.5 

0.64 ±0.08 

217 ±4 

85± 11 

48.0 ±3.3 

3900 ± 400 

0.15 ±0.02 

0.17 ±0.02 

Central 

1.0 

58.6 ± 7.0 

2.49 ±0.28 

236 ±2 

191 ±5 

42.5 ±2.0 

4500 ± 300 

0.58 ±0.05 

0.63 ±0.06 

West 

0.9 

22.8 ±2.7 

1.13 ±0.15 

258 ± 1 

113 ±8 

49.5 ±4.1 

3700 ± 400 

0.25 ±0.03 

0.29 ±0.02 


“ The root mean square noise levels were measured on either side of the H40a line at rest frequencies (based on v=236 km s“^) of 98.880-98.920 
and 99.135-99.175 GHz (excluding the unideintihed spectral feature at 98.915 GHz in the east region). 

^ These numbers are based on multiplying the measured continuum by 0.70 ± 0.10 to account for emission from sources other than free-free 
emission. 



Figure 3. The spectral energy d istribution of the cen t ral 30 arcsec of 
NGC 253 based on the data from lRodriguez-Rico et al.l (l2006h . The dark 
grey lines shows the function fit to the data, which comprises a synchrotron 
function with a power law index that was allowed to vary (shown in red), 
a free-free function that scales as (shown in green), and a modified 

Rayleigh-Jeans function for dust that varies as (shown in blue). 


recombination line emission using 
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which are adapted from IScoville & Murchikoval (l2013l) . The as 
term, which is the effec tive recombination coefficient listed by 
IStorev & Hummeii (Il995h . varies by a factor of ~3 between 5000 
and 15000 K but negligibly with electron density within the range 
10^-10^ cm“^. The SFRs listed in Table [T] are based on as cal¬ 
culated using Tie — 10^ cm“^. The two SFRs, which are mathe¬ 
matically linked through Equation [TJ are limited by both the pho¬ 
tometric accuracy and the assumptions behind the conversion be¬ 
tween Q and SFR. Because we are adjusting Te based on the line- 
to-continuum ratio, the calculated SFR from the continua vary by 
<lcr if we do not adjust the continuum emission to account for 
emission sources other than free-free emission. 

A comparison of these results to SFRs from ultraviolet, opti¬ 
cal, or near-infrared data would not be worthwhile, as the nuclear 
region is heavily obscured at those wavelengths (as discussed be¬ 
low). The mid- and far-infrared continuum emission tends to be 
either poorly resolved or saturated in most existing data, making 


it impossible to calculate star formation rates from dust emission 
on spatial scales comparable to the ALMA data. However, we can 
compare our total SFR (and Q, which is (3.2 zb 0.2) x 10^^ s“^ 
within the central 20x 10 arcsec) to equivalent measurements that 
have been made using millimetre and radio data (after rescaling all 
SFR and Q to correspon d to distances of 3.44 Mpc). 

IPuxlev et"^ (1 19971) obtained a distance-adjusted Q of (7.0 zb 
1.5) X 10^^ s“^ from H40a data that is significantly higher than 
our Q or other measurements of Q, but the line is detected at 
the ~5cr level, which may indicate i ssue s with the rel i ability of 
their data. iRodrfguez-Rico et al.l(l2006l) and lKenlev et al.l (l201ll) re¬ 
port Q based on lower-frequency recombination line data that are 
~3x lower than our measurement, his could be related to sensi¬ 
tivity issues since they had difficulty detecting line emission from 
all three sources that we detected. However, these groups may be 
probing primarily lower-density gas in the lower-frequency data, 
which would lead to lower estima tes of Q than what we obtain with 
H40a. lRampadarath et al.l ([201^ . in the latest 2.3 GHz analysis of 
supernova remnants in NGC 253, measure a SFR upper li mit of 
4.9 M© yr“^; our SFR falls below this limit. lOtt et al.l J200^ mea¬ 
sured a distance-adjusted SFR of 4.9 zb 0.5 M© yr“^ from 25 GHz 
continuum data which is ~3x higher than our value, but they rely 
upon a formula originating from ICondonI (Il992l) that uses a conver- 
sion of synchrotron emissio n to SFR may need to be recalibrated. 
iRodrfguez-Rico et al] (l2006h also produced an estimate of Q based 
on the continuum emission that is ~3x lower than our value, but 
as we indicated above, they may have biased the free-free emission 
estimates to low values. 

The optical recombination li ne emission from the cen tre of 
NGC 253 is heavily obscured, but lEngelbracht et ^ (Il998h mea¬ 
sured multiple near-infrared recombination lines, including Pa/3 
and Br 7 , that can be compared to H40a emission to infer the dust 
attenuation to the nucleus. If we assume that the dust functions like 
an attenuating sheet, then the dust attenuation A a at a given in¬ 
frared wavelength A is related to the ratio of the infrared line flux 
/(A) to the H40a line flux /(if40a) by 


Ax = —1.086 In 


eH40a 

3027 /xm 

r /(A) 1 

1 

A 

[fiH40a)\ 


(4) 


The lEngelbracht et al.l (Il998h spectroscopy data were acquired 
within a 12x2.4 arcsec aperture that also contains ~70% of the 
H40a e mission in our data, so we compare the Pa/3 and Br 7 fluxes 
from the lEngelbracht et al.l (Il998h spectroscopy data (which were 
the lines detected at the highest signal-to-noise levels) to 0.7 x 
our total H40a flux. Assuming that Ue = 10^ cm“^ and using 
Te = 6200 zb 400 K from Table [T] we obtain Apap — 5.0 zb 0.2 
and Apr-f = 4.2 zb 0 .2. This is higher than the values of Aj = 
2.00 zb 0.36 and Ak = 0.87 zb 0.16 for an attenuation sheet as re- 
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ported bv lEngelbracht et all jl99Sh . The results are affected by the 
assumed ne, although the effects are on the order of 0.2 dex for a 
range of 10^-10^ cm“^. 

It is likely that the attenuation was underestimated by 
lEngelbracht et al.l (Il998h . In their Br 7 image, they do not detect the 
same central peak that we detect in H40a (although this may be at¬ 
tributable to a coarser spatial resolution in the Br 7 data), and they 
measure a lower velocity dispersion than we do, implying that the 
emis sion from the centr al regio n is heavily obscured in the infrared 
data. lEngelbracht et al.l (Il998h indicate that the dust attenuation 
may be better described by a clumpy medium, which could affect 
the attenuation results. While it is beyond the scope of this paper to 
explore a more complex treatment of the attenuating medium, the 
comparison of Ax values based on the attenuat ing sheet scenario 
still su ggests that the attenuation measured by lEngelbracht et al.l 
( Il998ll is too low by ~3 dex. 

The broader implications of these results are that near-infrared 
recombination line emission is heavily obscured in very dusty 
galaxies, including starbursts like NGC 253 as well as luminous 
and ultraluminous infrared galaxies (galaxies with infrared lumi¬ 
nosities >10^^ L 0 ). SFRs based on near-infrared data may be 
underestimated in such objects. Near-infrared recombination line 
emission has been used to calibrate some SFR calculations based 
on combining ultraviolet/optical star formation tracers w ith mid- 
infrared star formation tracers (e.g. ICalzetti et al.l l2007h . While 
these formulae are probably still accurate for many galaxies, these 
equations may need to be recalibrated for application to more dusty 
objects. 


4 SUMMARY AND CONCLUSIONS 

Using ALMA observations of free-free and H40a emission at 
99.02 GHz from the centre of NGC 253, we obtain the following 
key results: 

• We measure Te of 3700-4500 K within the detected regions 
in the centre of the galaxy, which matches previous measurements 
within NGC 253 as is consistent with the range of values seen in 
the inner 8 kpc of the Milky Way. 

• Using both the continuum and line emission, we measure a 
Qof (3.2±0.2)xl0®'^ s“^ andaSFRofl.73±O.12M0yr“^ for 
the central 20x 10 arcsec. This measurement falls within the range 
of previously-published millimetre and radio results, which are in¬ 
consistent with each other, although our recombination line emis¬ 
sion is detected at superior sensitivity levels, and our constraints on 
the relative contribution of free-free emission at millimetre wave¬ 
lengths are more reliable. 

• The ratio of H40a emission to the near-infrared line emis¬ 
sion reported by lEngelbracht et al.l (Il998h implies that the central 
region is obscured by ~3 dex more than originally inferred. Given 
this, near-infrared line emission should be used with caution to 
measure SFR in very dusty regions or to calibrate other SFR met¬ 
rics. 

This analysis represents a sample of what can potentially be 
achieved using ALMA to measure Te and SFR. While this analysis 
is based on a target where millimetre and radio recombination line 
emission has been detected before, it is also based on data from 
ALMA Cycle 0 in which approximately one-third of the ALMA 
antennas were operational. Data from dedicated observations using 
ALMA when it is fully operational can potentially uncover mil¬ 
limetre recombination line emission in more galaxies, allowing us 
to probe star formation in sources where other data are not usable 
or to calibrate other SFR metrics. 
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